
Detailed Multifrequency Study of a Rapid VHE Flare of MrkSOl in May 2009 
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Abstract: We present observations of the gamma-ray blazar Markarian 501 between April 17 and May 5, 2009, with 
the Whipple 10-m telescope, VERITAS, the Fermi Large Area Telescope (LAT), Swift and RXTE as part of a 4.5-month 
multi-wavelength campaign. The presentation will focus on the strong very high energy (VHE) gamma-ray activity 
detected on May 1st with Whipple and VERITAS, when the measured flux {E > 400 GeV) reached five times the flux of 
the Crab Nebula, coincident with an increase in the optical polarization by a factor of 5, and a rotation of the polarization 
angle by 15 degrees. We also show that, during this 3- week period, the largest flux and spectral variability is seen at the 
highest energies of the broadband spectral energy distribution. 



Keywords: gamma rays, blazars, Whipple Telescope 

1 Introduction 

Blazars are a subclass of active galactic nuclei (AGN) with 
their relativistic jet pointing close to the line of sight. These 
objects are among the most powerful astrophysical sources 
known. They exhibit rapid and irregular variability, with 
extreme outbursts from minutes to days in all energy bands. 
The fast events at high energy are thought to be produced 
in a compact internal zone close to the system (e.g. lUl). 

For most blazars, the spectral energy distribution (SED) 
in a uFjy representation shows a double-peaked structure, 
where the first peak, from radio to X-ray, is attributed to 
synchrotron radiation and the second peak, in the GeV to 
TeV, to inverse-Compton or hadronic processes. 

Mrk 501 is a nearby 7-ray blazar (z = 0.034) first detected 
at TeV energies with the Whipple Telescope in 1996 12]. 
Since then, several multiwavelength campaigns were un- 
dertaken, mostly during episodes of high activity at long 
timescales or when flares on short timescales show up in 
the X-ray or 7-ray bands (e.g. [T], In 1997, the source 
presented the most intense episode of high activity detected 
at very high energies by the Whipple telescope reaching a 
level of 10 times the Crab Nebula flux, and it was observed 
in X-rays by BeppoSAX (L5J, 16J). 

2 Observations 

We present here observations of Mrk 501 from April 17 
to May 5, 2009, taken as part of a large-scale multi- 
wavelength campaign |7|, which included a number of 
ground and space-based experiments covering the spec- 



trum from radio to very high energy (VHE; £^ > 100 GeV) 
7-ray s. 

In this 3-week time interval, Mrk 501 was observed at the 
VHE with the Whipple Telescope [SJ every night for a total 
of 20 hours and with VERITAS ||9l[l0| for a total of 4 hours 
ifTHl . The VERITAS observations during the nights of April 
30 and May 1 were made with two telescopes, those on the 
remaining nights with three telescopes, because the other 
units of the full four-telescope array were affected by hard- 
ware issues. At the high energy (HE; E > 0.1 GeV) band, 
Fermi-hAI 1 12] monitored the source constantly, while the 
X-ray band was covered by Swift-XRT and RXTE-PCA. 
In optical and radio frequencies the source was covered by 
various instruments. In these proceedings we only show 
data from GASP, Mitsume and Swift-\J\OT at optical fre- 
quencies, and Metsahovi and OVRO at radio frequencies. 
More details of the instruments involved in the campaign 
and of the data analysis will be presented elsewhere. 

3 Results and discussion 
3.1 Light curves 

Figure 1 shows the light curves for some of the instruments 
involved in the campaign. In radio and optical bands, the 
measured fluxes are constant (within statistical errors). In 
the X-ray bands one can see statistically significant vari- 
ability (due to the high sensitivity of Swift-XKT and RXTE- 
PCA), but the flux variations are relatively small (well be- 
low 50 %). On the other hand, in the VHE domain, VER- 
ITAS and especially Whipple measured (statistically sig- 
nificant) flux variations of a factor of a few during this 
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Figure 1: Light curves for Mrk 501 from April 17 to May 5, 2009. Each dotted horizontal line represents a constant line 
fit for each instrument involved. Left: Top: OVRO at 15 GHz (filled squares) and Metsahovi at 37 GHz (open circles); 
Second: Swift-UYOT in the Ultraviolet, with three different bands, UVWl (260nm, open circles), UVM2 (220nm, filled 
stars) and UVW2 (193nm, open squares); Third: Mitsume in g band (nightly average); Bottom: GASP in the R band 
(nightly average). Right: Top: X-ray: Swift-XRT 0.3-2 keV (nightly average); Second: X-ray: RXTE-PCA (filled 
diamonds) and Swift-XRT (crosses) 2-10 keV (nightly average); Third: HE gamma-ray: Fermi-LKT (E > 300 MeV; 5- 
day average); Bottom: VHE 7-rays: Whipple {E > 400 GeV, normalised to E >300 GeV according to a power law with 
photon index -2.5; filled circles) and VERITAS (E > 300 GeV; filled squares; nightly average). 



time period and in particular up to a factor of 10 during 
MJD 54952. Therefore, during the above mentioned 3- 
week time interval, we find the highest variability at the 
highest energies. 

At VHE, the light curve is consistent with the constant 
emission of the source (3.9xl0~^^ph cm~^ s~^) until the 
night of May 1 (MJD 54952), when a high-emission state 
was detected with Whipple and VERITAS, reaching a max- 
imum 7-ray flux 10 times the average baseline flux, which 
is approximately 5 times the Crab Nebula flux. 

Figure 2 shows the Whipple 10-m and VERITAS fight 
curves for May 1, 2009 where each point corresponds to 
a 4-minute bin. The flux increased by a factor of ~ 5 in the 
first 30 minutes. In the days after the flare (MJD 54953- 
55), the source continued in a high state, and the flux each 
night was about twice the baseline flux. 

The Swift-XKT observations are not strictly simultaneous 
to the VHE observations performed by Whipple and VER- 
ITAS. The time difference between Swift observations and 
Whipple and VERITAS observations is 7 hours. Given 
that this source is known to have very fast variability (e.g., 
see f3l), not having strictly simultaneous observations is a 
caveat that cannot be neglected. In any case, the VERI- 
TAS observations started 1.5 hours after the Whipple ob- 
servations and continued with simultaneous (to Whipple) 
observations until the end of the night. These two instru- 
ments observed a flux enhancement by more than a factor 
of 5. Therefore, it is reasonable to assume that the VHE 
flare lasted more than 3 hours, and hence that Swift's obser- 
vations might have occurred during this bright VHE state. 
However, Swift did not record any substantial flux increase 



during MJD 54952. The flux in the energy range 0.3-2 
keV is essentially compatible with the flux during previ- 
ous days. In the range 2-10 keV there is a flux increase on 
MJD 54952, which lasted several days. However, the mag- 
nitude of the flux increase in the 2-10 keV band is about 
50%, which is substantially smaller than the flux increase 
observed by Whipple and VERITAS at the VHE domain. 

Fermi-LAI operates in a survey mode, which each point 
of the sky is observed approximately during 30 min every 
3 hours and hence sources are observed continuously on 
timescales down to 3 hours. However, Mrk 501 is a rela- 
tively weak source for Fermi-LAI and hence one typically 
needs to integrate over several days in order to have a sig- 
nificant detection. In these proceedings we present a light 
curve with the data binned in 5 -day time intervals. The 
second-to-last time interval starts on MJD 54952 and con- 
tains the entire VHE flare, but we do not see any signifi- 
cant variation with respect to the previous time intervals. 
However, it is worth mentioning that on MJD 54952 (day 
with the highest VHE flux), Fermi-LAI detects Mrk 501 in 
a single day with a TS larger than 25 (i.e. a significance 
greater than ~ 4. 6a), which does not occur for the other 
days. The measured flux (above 300 MeV) for this day is 
(3.46 ± 2.37) 10-^ ph cm-^ 

Optical observations were performed as part of the Steward 
Observatory balzar monitoring program with the (2.3 m) 
Bok and the (1.54 m) Kuiper telescopes every night from 
MJD 54947 to MJD 54955, which includes the night of the 
VHE flare. Figure 3 shows the light curve, the degree of the 
optical linear polarization and the electric-vector position 
angle (EVPA) from those nights. 
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Figure 2: Whipple 10-m and VERITAS light curve (4- 
minute binning) for the night of the VHE flare. The solid 
line shows the baseline emission for the source depicted in 
Figure 1. 
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Figure 3: Optical flux, degree of the optical linear polar- 
ization and EVPA light curves (first, second and third plots 
respectively) measured with the Steward Observatory, and 
(forth plot) the VHE light curve obtained with Whipple and 
VERITAS. 



From MJD 54947-51, the degree of polarization was 
steady, and dropped a 15 % after the VHE flare. The 
EVR\ light curve shows a continuous increase from 15 to 
30 degrees in 3 days, and the rotation stops right when the 
large VHE flare occurs. If the two events can be physically 
linked, that could indicate a common origin for the optical 
and 7-ray emission, as already seen in |[T3l[T4l[T5]i 

3.2 Variability & Correlation 

In order to quantify the flux variability in the light curves, 
the fractional RMS variability amplitude, Fyar fl^l was 
calculated as: 
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Figure 4: Fractional variability amplitude for all the instru- 
ments involved. 



(1) 



where F is the average photon flux, S the standard 
deviation of the N flux measurements and is the mean 
squared error. The uncertainty of Fyar is given by: 



AF,a 



A 




(2) 



These values were calculated using a daily average for each 
energy band. Figure 4 shows the Fyar values obtained for 
all the experiments involved, except for Fermi-LAI where 
the result showed a negative excess variance (a^ > S'^), 
indicating low level of variability and/or slightly overesti- 
mated errors. Essentially such a result can be interpreted 
as no signature for variability either because there was no 
variability or the instrument was not sensitive enough to 
detect it. 

The value of Fyar is compatible with zero or very low for 
all the energy bands, except for the VHE range, where 
VERITAS observed a Fyar of 0.62±0.05 and Whipple 
0.95 ±0.07. This large variability in the VHE domain is 
clearly dominated by the large VHE flare observed on MJD 
54952 and the few days after. 

Historically, the X-ray and VHE 7-ray fluxes are correlated 
in general, yet we found that the large VHE 7-ray flare was 
not accompanied by a large X-ray flare. The Swift-XKT ob- 
servations are not simultaneous, but are contemporaneous 
(7 hours late) and did not show any increase in the flux. 
The flux in the 2-10 keV band showed an increase of about 
50% during the days following the VHE flare, but this flux 
enhancement is substantially smaller than the one seen by 
Whipple and VERITAS during the same days. At the GeV 
band the source remained steady with no variations during 
those days. To investigate the relationship between the X- 
ray and the VHE 7-ray, the discrete correlation function 
(DCF) was calculated (Tl] and does not indicate a signifi- 
cant correlation between the two datasets. We can conclude 
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Table 1 : Best-fit parameters for VHE data. 
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Figure 5: Time-averaged VERITAS and Whipple photon 
spectra of Mrk 501 for discrete flux levels (see text). 



that this VHE flare is very probably an "orphan" flare (in 
the context of not being accompanied by an X-ray flare), 
like the one previously detected on lES 1959+650 118 |. Al- 
ternatively, the X-ray flare could exist, but being substan- 
tially smaller than that seen at VHE, as it occurred for PKS 
2155-304 in 2006 1 19|. Further details on this subject wiU 
be given on a forthcoming publication. 

3.3 Spectral Energy Distribution 

The evolution of the TeV energy spectrum during the short 
flaring state of Mrk 501 could be important for under- 
standing the mechanism of particle acceleration in a blazar 
source. The differential energy spectra of Mrk 501 (see 
Figure 5) with the Whipple Telescope and VERITAS were 
modeled for the quiescent emission and for the flare state 
with a simple power law in each case: 



=FoxlO-^(£'/lTeV)-rvHE m- 



TeV- 



where Fq is a normalization factor and Fvhe is the photon 
index. The best-fit parameters and associated errors for the 
VHE data are summarized in Table [T] 

A marginal indication of spectral hardening with increas- 
ing flux activity is found in the TeV band. A similar trend 
had already been found in earlier observations in 2005 with 
MAGIC |3 1, when the softest photon index was 2.43±0.05 
(for the low/medium state), and the hardest photon index 
was 2.09±0.03 (for the flare on MJD 53551). 

The broadband spectral results obtained for the high state 
(MJD 54952-55) of Mrk 501 during the three-week period 
(MJD 54936-56) are plotted in Figure 6. For comparison, 
the average spectrum for the 4. 5 -month multi wavelength 
campaign (which includes the three-week period) is also 
plotted |7|. Figure 6 shows very clearly that the largest 
variation in the emitted flux occurs in the VHE domain. 
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Figure 6: Spectral energy distribution of Mrk501 during the 
interval MJD 54952-54955 (depicted with red circles), and 
the VHE spectrum measured by Whipple on MJD 54952 
(depicted with green stars). For comparison purposes, the 
4. 5 -month averaged SED from the entire campaign (ex- 
tracted from 1 7 1) is also shown with black diamonds. 
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